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Abstract 

This paper reports the discovery of a bright X-ray transient source, Suzaku J1305- 
4913, in the south-west arm of the nearby Seyfert II galaxy NGC 4945. It was detected 
at a 0.5 - 10 keV flux of 2.2 x 10~^^ erg cm~^ s~^ during the Suzaku observation 
conducted on 2006 January 15 - 17, but was undetectable in a shorter observation on 
2005 August 22 -23, with an upper limit of 1.7 x 10"^^ erg cm~^ s~^ (90% confidence 
level). At a distance of 3.7 Mpc, the bolometric luminosity of the source becomes 
Lboi = 4.4 X lO^^a erg s~^, where a = (cos 60°/ cosz) and i is the disk inclination. 
Therefore, the source is classified into so-called ultraluminous X-ray sources (ULXs). 
The time-averaged X-ray spectrum of the source is described by a multi-color disk 
model, with the innermost accretion disk temperature of Tin = 1-691q;o5 keV. During 
the 2006 January observation, it varied by a factor of 2 in intensity, following a clear 
correlation of Lboi oc T^. It is inferred that the innermost disk radius i?in stayed 
constant at -Rin = 79tl'^a^^'^ km, suggesting the presence of a standard accretion disk. 
Relating i?in with the last stable orbit around a non-rotating black hole yields a 
rather low black hole mass, ~ 9a^/^ solar masses, which would imply that the source 
is shining at a luminosity of ~ 3a^^'^ times the Eddington limit. These results can be 



better interpreted by invoking sub-Eddington emission from a rapidly spinning black 
hole with a mass of 20 - 130 solar masses. 

Key words: galaxies: individual (NGC 4945), X-rays: galaxies, black hole 
physics, accretion, accretion disks 

1. Introduction 

A number of nearby galaxies host bright off-center point-like X-ray sources, with lumi- 
nosities of 10'^^ - 10^° erg s~^ (e.g., Fabbiano 1989). These X-ray sources are usually called 
ultraluminous X-ray sources (ULXs; Makishima et al. 2000). Based on their high luminosities, 
ULXs are regarded as candidates for intermediate-mass black holes with masses of M = 20 - 
1OOM0. X-ray spectra of ULXs are usually described by a multi-color disk (MCD; Mitsuda et 
al. 1984; Makishima et al. 1986) model, or a power-law (PL) model. Because the MCD model 
provides a good approximation for the X-ray emission from a standard accretion disk (Shakura 
& Sunyaev 1973) and successfully interprets the observed X-ray spectra of Galactic black hole 
candidates, the black hole scenario for ULXs seems plausible. However, a simple application 
of the MCD model to ULXs gives an innermost disk temperature (Tin = 1-2 keV, Makishima 
et al. 2000) significantly higher than those of Galactic black holes (Tin = 0.5 - 1 keV). This 
is apparently inconsistent with the argument for a intermediate-mass black hole, since a more 
massive object would have a lower disk temperature M According to recent stud- 

ies (e.g., Mizuno et al. 2001; Watarai et al. 2000; Mizuno et al. 2007; Tsunoda et al. 2006), 
the high temperature of ULXs can be explained if they harbor "slim disks" (Abramowicz et 
al. 1988; Szuskiewicz et al. 1996; Watarai et al. 2000), in which an advection-dominated accre- 
tion flow and/or photon trapping (e.g., Ohsuga et al. 2005) become important, due to a high 
accretion rate. 

Among accreting Galactic black holes that are accompanied by standard accretion disks, 
some objects exhibit rather high temperatures, such as GRO J1655-40 and GRS 1915+105. 
Zhang et al. (1997) proposed a spinning (i.e., Kerr) black hole scenario, for such objects; the 
last stable orbit of a prograde accretion disk around a Kerr black hole can get closer to the black 
hole. Accordingly Tin can potentially be higher, compared to that of a non-rotating black hole, 
although various relativistic effects should be taken into account. The Kerr black hole scenario 
was subsequently applied to ULXs by several authors (e.g., Mizuno et al. 1999; Makishima et 
al. 2000). Moreover, recent X-ray observations, including those with Suzaku in particular, have 
revealed evidence of black hole rotation in the Seyfert I galaxy MCG— 6— 30— 15 (e.g., Miniutti 
et al. 2007). 

NGC 4945 is a nearby spiral galaxy with an edge-on configuration (~ 80°; Ott et al. 
2001). It is regarded as a member galaxy of the Centaurus cluster, and its distance is precisely 
estimated to be 3.7 Mpc (Mauersberger et al. 1996). This galaxy hosts a Seyfert II nucleus 

2 



which is one of the brightest known active galactic nuclei in the universe above 20 keV (Done 
et al. 1996; Guainazzi et al. 2000; Done et al. 2003). While the nucleus is directly invisible due 
to heavy obscuration below 10 keV (Iwasawa et al. 1993), several bright X-ray point sources 
are detected within the galaxy (Guainazzi et al. 2000; Schurch et al. 2002), some of which are 
classified as ULXs (Feng & Kaaret 2005). 

We observed NGC 4945 twice with Suzaku (Mitsuda et al. 2006) separated by 5 months, 
and detected a bright transient ULX in its south-west arm region, during the second observation 
in January 2006. The observed spectral behaviour of the source indicates that it harbors a 
standard accretion disk rather than a slim disk, but simply invoking a Schwarzschild black 
hole leads to an implausible conclusion that a standard disk is shining at a significantly super- 
Eddington luminosity. Therefore, we have adopted the spinning black hole scenario, and have 
confirmed that the observed properties can be consistently explained as sub-Eddington emission 
from a rapidly rotating black hole, with a mass of 20 - 130 Mq. 

2. Observation and Data Reduction 

The present observation of NGC 4945 with the fifth Japanese X-ray observatory Suzaku 
(Mitsuda et al. 2006) were performed twice, during the scientific working group (SWG) phase; 
the first was conducted on 2005 August 22 - 23 for about 27 ks, and the second one on 2006 
January 15 - 17 for about 100 ks. The X-ray Imaging Spectrometer (XIS; Koyama et al. 2006) 
and the Hard X-ray Detector (HXD; Takahashi et al. 2006; Kokubun et al. 2006) onboard 
Suzaku were operated in the normal clocking mode with no window option and in the normal 
mode, respectively. The nucleus of NGC 4945 was placed at the XlS-nominal and HXD-nominal 
positions (Serlemitsos et al. 2006) in the first and second observations, respectively. 

In the present paper, we concentrate on the XIS data; the HXD data on the nucleus are 
presented separately by Itoh et al. (2007). We analyzed the data with Revision 0.7 processing 
(an internal release for the SWG; Fujimoto, et al. 2006), with the HEADAS 6.0.6 software 
package. Following standard data processing, we filtered the data by adopting the criteria that 
the spacecraft is outside the south Atlantic anomaly, the geometric cut off rigidity is higher 
than 6 GV, the source elevation above the rim of bright and night Earth is higher than 20° 
and 5° respectively, and the XIS data are free from telemetry saturation. As a result, about 
23 and 80 ksec of good exposures have remained from the first and the second observations, 
respectively. In the scientific analysis below, we utilize only the events with a grade of 0, 2, 3, 
4 or 6. 

Because we noticed some systematic uncertainties in the position determination in these 
observations, we corrected the Suzaku attitude solution, utilizing the archival XMM-Newton 
data (Obsid: 0204870101). By referring to the positions of X-ray sources detected with XMM- 
Newton (crosses in Figure 1) in a region outside NGC 4945, but within the XIS field of view, 
we shifted the coordinates of the XIS data by about 15" and 45", for the first and second 
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observations, respectively. These angular offsets are within nominal attitude tolerance of the 
Revision 0.7 of Suzaku data (about 1'; Serlemitsos et al. 2006). 

3. X-ray Image 

Figure 1 shows 0.5 - 10 keV X-ray images of NGC 4945, obtained with the XIS in the 
two observations. Only data from the three front-illuminated (FI) CCD sensors (XISO, 2, and 
3; Koyama et al. 2006) are summed up. An infra-red (IR) image of NGC 4945 (Jarrett et al. 
2003) is superposed with contours in both observations. In addition to some point sources, we 
detected diffuse X-ray emission, elongated along the IR image of the galaxy. The X-ray peak 
of NGC 4945 coincides with the nucleus of NGC 4945, within 8". 

Although the two Suzaku images are similar, a clear difference is a bright transient X- 
ray source which has emerged in the second observation in the south-west arm of NGC 4945. 
After the attitude correction, the X-ray peak position of the transient source was determined 
to be {a, 6) = (13^05™05f5, -49°31'39") in J2000 coordinates, with an accuracy of about 10". 
The source is found to have no counterpart among the XMM-Newton sources shown in figure 
1 (crosses), of or recent X-ray catalogs (Ptak & Griffiths 2003; Colvert et al. 2004). With 
ROSAT, Brandt et al. (1996) reported the detection of an off-center X-ray transient source, 
NGC 4945 X-1. However, the position of NGC 4945 X-1, {a, 6) = (13^05"'22f6, -49°29'12") 
with an error of 15", is clearly different from that of our Suzaku source. The two transients 
are thus considered to be distinct. Therefore, we regarded this Suzaku source as a new X-ray 
source, and named it Suzaku J1305-4931. 

In 2005 February 8, a supernova, SN 2005aj, was optically discovered (Jacques & 
Pimentel 2005) in the south-west of NGC 4945 near Suzaku J1305-4931. However SN 2005af, 
[a, 6) = (13^'04'"44?06,-49°33'59:'8), was outside the XIS field of view in both observations. 
Therefore, Suzaku J1305-4931 cannot be the X-ray counterpart of this supernova, either. In 
the following, we assume Suzaku J1305-4931 to be physically associated with NGC 4945. 

4. The X-ray properties of Suzaku J 1305-4931. 

4.I. Set up for analysis 

In order to investigate the properties of Suzaku J1305-4931, we integrated source and 
background X-ray signals within circular regions, which are denoted as SRC and BGD in the 
right panel of figure 1, respectively. Both regions have a radius of 1.'67, or 1.79 kpc in physical 
size at the distance of 3.7 Mpc. It is difficult to adopt a larger integration circle, because Suzaku 
J1305-4931 is detected near the edge of the XIS field of view. The BGD region has the same 
separation from the NGC 4945 nucleus as the SRC region. 

In analyzing the XIS spectra, we adopted a response matrix function 
ae_xi [0,2,3] _20060213.rmf for the FI CCD sensors, and ae_xil_20060213c . rmf for 
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the backside-illuminated (BI) sensor (XISl; Koyama et al. 2006), all of which are appropriate 
for the revision 0.7 data. We discarded the 1.82 - 1.84 keV region, where the XIS response is 
not yet fully calibrated due to the presence of Si-K edge. 

Due to vignetting effects of the X-ray telescope feeding the XIS, the effective area for the 
SRC region in the second observation is ~ 40 % smaller than that for the XIS or HXD nominal 
position. In addition, the SRC region (1.'67 in radius) is significantly smaller then that of typical 
integration regions for a point source (e.g. ~ 3'), so that the effective area is further reduced. 
Therefore, we calculated the auxiliary response files for both observations using a Monte-Carlo 
simulator xissimarfgen (Ishisaki et al. 2006), as of 2006 May 26 (an internal release for the 
SWG). Soft X-ray absorption due to the contaminants on the optical blocking filters over the 
XIS CCD cameras (Koyama et al. 2006) were taken into account in the simulation. 

4-2. X-ray lightcurves 

Figure 2 shows XIS FI lightcurves of Suzaku J1305-4931 in the soft (0.5 - 2 keV) and 
hard (2 - 10 keV) energy bands, accumulated over the SRC region. The background has been 
subtracted utilizing the BCD region. After summing over the three XIS FIs, we measured the 
time-averaged net count rate of the source in the soft and hard bands to be (8.3 ±0.1) x 10"^ 
cts s~^ and (8.7 ±0.1) x 10"^ cts s~^, respectively. The source thus varied by a factor of 2 
on a time scale of about half a day in both energy bands, with the hard band exhibiting a 
larger amplitude. The bottom panel of figure 2 shows the hardness ratio, simply calculated by 
dividing the hard band count rate by the soft band one. It suggests that the hardness of the 
source changed with its brightness. We further examine the spectral variation in §4.5. 

4-3. X-ray emission from the host galaxy. 

Before analyzing the X-ray spectrum of Suzaku J 1305-4931, we briefly estimate the 
contribution from the NGC 4945 galaxy itself to the SRC region. For this purpose, we extracted 
the XIS spectra from the first observation, and show them in figure 3. The data were binned 
into energy intervals each containing at least 30 events. 

We fitted the XIS spectra with a two component model, consisting of a soft thermal 
emission (a MEKAL model; Mewe et al. 1985; Mewe et al. 1986; Kaastra et al. 1992; Liedahl et 
al. 1995) and a PL model. The former and latter are meant to represent emission from thermal 
plasmas associated with the galaxy, and the integrated emission from unresolved point sources 
within the galaxy (including two faint XMM-Newton sources; figure 1), respectively. Both 
components were subjected to the Galactic absorption toward NGC 4945, with a hydrogen 
column density of A^h = 1-6 x 10^^ cm~^ (Dickey & Lockman 1990). As summarised in Table 
1, we obtained an acceptable fit with reasonable parameters (e.g., Guainazzi et al. 2000). 
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4-. 4- Time-averaged X-ray spectrum 

Figure 4 shows the XIS spectra of Suzaku J1305-4931, obtained by accumulating events 
from the 2006 data over the SRC region, and then subtracting background spectra from the 
BGD region. The FI and BI spectra were binned into energy intervals with at least 100 and 
50 events, respectively. While the X-ray signals are significantly detected in 0.5 - 10 keV, we 
discarded the data below 0.75 keV in the spectral fitting, due to a severe contribution by the 
emission from the host galaxy (figure 3). The spectrum appears relatively featureless, without 
any emission nor absorption lines. We fitted the spectrum first with a PL model, and then with 
an MCD model; both models are widely adopted in studying X-ray spectra of Galactic black 
hole binaries and ULXs. In the fitting, we took into account the X-ray emission from the host 
galaxy by the best-fit two-component model for the X-ray spectrum of the SRC region in the 
first observation (Table 1). Since we noticed a calibration uncertainty in the flux determination 
between XIS FI and BI, we allowed in the fittings the relative model normalization to deffer 
between FI and BI. As a result, XIS BI was found to yield a flux lower by about 15 % than 
XIS FI. Throughout the present paper, we refer to the XIS FI flux. 

The PL model gave a best-fit photon index of F = 1.89lo;o5, with a moderate fit goodness 
(x^/d.o.f. = 278.8/233). However, as shown in the second panel of figure 4, the fit residuals 
(particularly above 2 keV) imply that the observed spectra have a somewhat more convex shape 
than a single power law. This motivated us to attempt an MCD model fit. 

The MCD model gave a slightly better fit (the third panel of figure 4, x^/d.o.f. = 
273.8/233) than the PL one, although the MCD model appears to be slightly too convex. 
The absorption column density, A^'h = 2.17to;26 ^ 10^^ cm~^, is larger by about 30 %, than 
the Galactic value toward NGC 4945. The innermost temperature and innermost disk radius 
were measured to be Tin = l-69lo;o5 keV and Ri^ = 77.0l|;7 a^/^ km, respectively, where a = 
(cos 60°/ cosz) with i being the inclination of the accretion disk to our line of sight. Here, we 
converted the apparent radius Tin to the true radius as i?in = K^^rin, by adopting two correction 
factors (Makishima et al. 2000 for summary); a spectral hardening factor of k ~ 1.7 (Shimura & 
Takahara 1995), which represents the ratio of the color temperature to the effective temperature, 
and a correction factor for the inner boundary condition, ^ = 0.412 (Kubota et al. 1998). The 
bolometric luminosity of the accretion disk was calculated to be Lboi = 4.4 x lO^^a erg s~^, using 
the equation Lboi = 47r(i?in/,^)^(7(Tin/ k)^, where a is the Stefan-Boltzmann constant. This value 
of Lboi classifies Suzaku J1305-4931 into the ULX category. Figure 5 shows the spectral energy 
distribution of Suzaku J1305-4931, together with the best-fit MCD model, in the uF,^ form. 

We also fitted the observed spectrum with the MCD + PL model. As shown in figure 6 
and Table 2, the fit was significantly improved (x^/d.o.f. = 235.2/233) by this composite model. 
The PL component dominates the MCD one above ~ 5 keV, and compensates the spectral 
curvature more concave than the sing le MCD model. The MCD temperature, = 1.25l|5:?f 
keV, becomes slightly lower than that of the single MCD fit. However, the PL component 
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turned out to be rather hard, T = 1.23_i;47, compared with the typical value of ULXs with 
MCD+PL-type X-ray spectra (F > 2.0; e.g, Colbert & Mushotzky 1999; Winter et al. 2006). 
Therefore, we regard this composite model as physically of little meaning, and do not discuss 
it in any further. Incidentally, this MCD + PL model fit did not converge to a "cool" disk 
(e.g.. Tin ~0.2 keV) disk solution, which was reported by some authors (e.g.. Miller et al. 2003). 
Likewise, based on the Suzaku data on NGC 1313 X-1, Mizuno et al. (2007) cast a doubt upon 
the reality of such "cool" disks in ULXs. 

4-5. X-ray Spectral Variability 

Since figure 2 suggests the intensity-correlated spectral variation of Suzaku J 1305-4931, 
we plotted in figure 7, the hardness ratio against the total (0.5 - 10 keV) count rate. The 
plot shows a clear trend that the spectrum of the source gets harder, as it becomes brighter. 
Therefore, we divided the total exposure into 4 periods according to the intensity of the source, 
and investigated in detail the variation of the spectral parameters. The definition of these 
periods are indicated with dashed- lines in figure 2, and their effective exposures are shown in 
table 3. 

Figure 8 shows the XIS spectra of these 4 periods. We fitted them by the MCD model, 
because it gave a better fit to the time-averaged spectra than the PL model. The absorbing 
column density was fixed at the best-fit value from the time-averaged spectra shown in table 
2, A^H = 2.17 X 10^^ cm~^. The model successfully described the four spectra with the best-fit 
parameters tabulated in table 3. The source shows the highest temperature as Tin = l-77lo;o8 
keV (i.e. the hardest spectrum) in the brightest period (Period 1, Tboi = 5.9 x lO'^^a erg 
s~^), while the lowest temperature as Tin = l-Sltan keV in the faintest period (Period 4, 
Tboi = 3.0 X lO'^^a erg s~^). This is consistent with the behaviour of the hardness ratio shown 
in figures 2 and 7. 

Throughout the observation, the innermost radius of the accretion disk has been found 
to stay relatively constant at ~ SOa^/^ km. Therefore, we re-fitted the spectra of the four 
periods, by tying R^^ together, but allowing Tin to take independent values. We also adopted a 
common relative normalization factor between XIS FI and BI. A similar fit goodness has been 
obtained (x^/d.o.f. =485.0/516), in comparison with a simple sum over the four individual fits 
(x^/d.o.f. = 484.4/510 in total); the other parameters remained unchanged within statistical 
uncertainties. As a result, the disk radius has been obtained as i?in = 79.lj;3;g a^^^ km, which 
agrees, within errors, with that from the time-averaged spectrum. 

5. Discussion 

5.1. Summary of the results 

In the second Suzaku observation of the nearby Seyfert II galaxy NGC 4945, performed 
on 2006 January 15 - 17, we have discovered a new bright transient ULX, Suzaku J1305-4931, 
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within the south-west arm region. We assume the source to be physically associated with NGC 
4945. The time- averaged XIS spectrum of the source was approximated by a PL model with 
photon index ~ 1.9, or somewhat better by an MOD model with Tin ~ 1-7 keV. Adopting the 
MCD model, the bolometric luminosity of the source was determined to be Lboi = 4.4 x lO^^a 
erg s~^. These values lie within the typical range of those of well-studied ULXs (e.g. Makishima 
et al. 2000; Mizuno et al. 2001). During the observation, Suzaku J1305-4931 varied by a factor 
of 2 in intensity in a timescale of about a half day. 

Suzaku J1305-4931 was undetectable in the first Suzaku observation, conducted in 2005 
August 22 - 23, and has no counterpart among the available X-ray source catalogs (Ptak & 
Griffiths 2003; Colvert et al. 2004). We fitted the XIS spectrum from the first observation, by 
adding an MCD component to the emission from the host galaxy (Figure 3 and Table 1). Then, 
the upper limit on the 0.5 - 10 keV source fiux was obtained to be 1.7 x 10"^^ erg cm~^ s~^ (at 
90% confidence level). Therefore, the source is inferred to have brightened by more than two 
orders of magnitude in about 5 mounts. 

The MCD model fit gives an innermost disk radius of i?in = 77.0^4^7 a^^^ km. Equating 
i?in to the last stable orbit of a standard accretion disk (Shakura & Sunyaev 1973) around a 
Schwarzschild black hole, SRg where Rs is the Schwarzschild radius, the mass of the source 
would become M = (i?in/8.86 km)M0 = 8.7toia^^'^MQ. However, the Eddington limit for this 
mass is calculated as i^Edd = 1-3 x lO'^^a^/^ erg s~^. Therefore, the source would be concluded 
to be shining at the Eddington ratio oi rj = Lboi/-^^Edd = source appears to 

show the "super-Eddington" problem. 

5.2. Slim disk or standard disk ? 

The "super-Eddington" problem is commonly found among ULXs showing an MCD type 
X-ray spectrum. Usually, the problem is solved by invoking the slim disk model (e.g, Mizuno 
et al. 2001; Watarai et al. 2000; Tsunoda et al. 2006), which is thought to be realized at a 
high accretion rate (Abramowicz et al. 1988; Szuskiewicz et al. 1996; Watarai et al. 2000). In a 
slim disk, an advection-dominated accretion fiow, and/or photon trapping (e.g., Ohsuga et al. 
2005) become important. As a result, the X-rays are expected to be radiated from inside the 
last stable orbit of the accretion disk as well, and the local temperature profile is predicted to 
become fiatter than the standard disk (Watarai et al. 2000). 

In order to examine the present data for the possible presence of a slim disk, we ana- 
lyzed the time-averaged XIS spectrum with an accretion disk model, in which the local disk 
temperature is proportional to r"^ with p being a positive free parameter, where r is a local 
radius (Mineshige et al. 1994; Hirano et al. 1995; Kubota & Makishima 2004). While the MCD 
model assumes the fixed value of p = 3/4, integrated emission from a slim disk is thought to be 
approximated by smaller value of p, down to p = 1/2 (Watarai et al. 2000). Here, we denote 
this model the "variable-p" disk model. A good fit was obtained by this model (the forth panel 
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of figure 4; x^/d.o.f. = 243.5/232), together with the parameters shown in table 2. Compared 
with the MCD fit, the disk temperature became higher. Tin = 2.3llo;23 keV, and the disk radius 
got smaller, _Rin = 37.6li3;o <^^^^ km. We obtained p = 0.59 ± 0.03. This value of p is consistent 
with the fact that the observed X-ray spectrum is slightly more convex than the PL model, 
but slightly more concave than the MCD one. All these parameters appear to be reasonable 
for the slim disk picture. 

The observed spectral variation of Suzaku J1305-4931, on the other hand, seems to 
contradict the slim disk interpretation. In Figure 9, we plot Lboi against Tin of the source, 
for the 4 individual periods (see Table 3). For comparison, we also compiled those of some 
well-studied Galactic and Magellanic black holes, and ULXs (Makishima et al. 2000; Kotoku 
et al. 2000; Mizuno et al. 2001; Miyawaki et al. 2004; Tanaka et al. 2005; Tsunoda et al. 2006). 
The super-Eddington problem of Suzaku J1305-4931 is clearly visualised in the figure, in which 
the data points for the source are located above the line of 77 = 1. However, it is important 
to note that Suzaku J1305-4931 shows a clear correlation as Tboi oc T^. Correspondingly, the 
innermost disk radius stays constant around -Rin = 79.ll3;9 a^^^ km, which is consistent with 
the value of the time-averaged spectrum. Such correlations have been frequently observed in 
Galactic/Magellanic black holes (e.g.,Ebisawa et al. 1993; Kubota & Makishima 2004), and is 
considered to be a signature of a standard accretion disk (Shakura & Sunyaev 1973). On the 
other hand, the MCD model fit to the numerically simulated slim disk spectra is reported to 
give a relation of Lboi ocT^^ (Watarai et al. 2000), and various ULXs actually follow this relation 
(Figure 9; Mizuno et al. 2001). 

For a further examination, we simulated the dependence of the hardness ratio on the 0.5 
- 10 keV count rate for various values of Tin, assuming either Lboi oc or Tboi oc T^. We utilized 
the same response files as used in the spectral fitting, and took into account the contribution 
of the emission from the host galaxy (table 1). The model predictions are shown in figure 7, in 
comparison with the observed data. The standard disk (Tboi oc T-^^) locus is found to describe 
the observed relation between the hardness and count rate with x^/d.o.f. = 9.9/11 (only the 
erros of the hardness ratio are considered), which is better than is obtained with the slim disk 
(Lboi oc T^) one x^/d.o.f. = 12.7/11. Therefore, the source variation is described better by the 
Lboi oc relation than by the Lboi oc T^^ dependence. 

Except for Galactic/Magellanic sources, a luminous source, NGC 253 X21 (Tanaka et 
al. 2005), with a luminosity of Tboi ^ 10'^^ erg s~^, is reported to follow a similar Tboi oc 
spectral behavior. However, this correlation has not been found so far from ULXs with a 
luminosity higher than a few 10^^ erg s~^. Therefore, we regard Suzaku J1305-4931 to be the 
most luminous black hole that shows the Lboi oc T-^ relation, and hence is thought to host a 
standard accretion disk, although a simple standard disk cannot solve the super-Eddington 
problem. 
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5.3. Kerr black hole scenario 

We have to look for some other mechanisms to explain the high luminosity and high 
temperature of the source. One of the most important remaining physical parameters is a 
possible rotation of the black hole. The rotating black-hole scenario (i.e., a Kerr black hole) 
was originally proposed by Zhang et al. (1997) to explain the high disk temperature of the 
Galactic jet source, GRS 1915+105 and GRO J1655-40. The idea was subsequently applied 
to the interpretation of the X-ray properties of the ULXs observed with ASCA (Mizuno et al. 
1999; Makishima et al. 2000). Rapid black-hole spin is expected to make the last stable orbit of 
an accretion disk around it come closer to the black hole, down to 0.5-Rs for the most extreme 
case. Neglecting the general relativistic effects, the Kerr black hole assumption can enlarge 
the mass estimated by equating R^^ to the last stable orbit, by up to 6 times. Simultaneously, 
the maximum value of Tin at the Eddington accretion rate can increase by a factor of -\/6- 
Moreover, the radiation efficiency of a Kerr black hole can be significantly higher than that 
of a Schwarzschild black hole (Zhang et al. 1997). Therefore, the Kerr black hole scenario is 
thought to give a possible solution to the super Eddington problem of Suzaku J1305-4931. 

In order to quantitatively evaluate the possibility of black hole rotation in Suzaku J1305- 
4931, we adopted the KERRBB model in the XSPEC package (Li et al. 2005). The model 
describes the numerically computed radiation spectrum from a geometrically thin accretion disk 
around a Kerr black hole, and adequately takes into account various general relativistic effects. 
It is intended to extend the GRAD model in XSPEC (Hanawa 1989; Ebisawa et al. 1991), 
which represents the general relativistic accretion disk around a Schwarzschild black hole. Li et 
al. (2005) applied this model to three Galactic sources, for which the mass and disk inclination 
are already known; GRO J1655-40, 4U 1543-47, and XTE J1550-564. They concluded that 
GRO J1655-40 and 4U 1543-47 have a high black-hole spin, while XTE J1550-564 hosts a 
Schwarzschild black hole; the results on GRO J1655-40 is consistent with a speculation by 
Kubota et al. (2001) based on its comparison with LMC X-3. Therefore, the KERRBB model 
is quite suitable for our purpose. 

We re-analysed the time-averaged spectrum of Suzaku J1305-4931 with KERRBB. 
Among the nine physical parameters in the model, we switched off the effects of limb-darkening, 
self-irradiation and torque at the inner boundary of the accretion disk, since these are relatively 
unimportant for our purpose (Li et al. 2005). We again adopted a spectral hardening factor 
of K = 1.7 (Shimura & Takahara 1995). Then, the emergent spectrum should be determined 
by specifying the mass M, mass accretion rate M, disk inclination i, and spin parameter 
a = c/{GM'^)J, where c, G, and J are the light speed, the gravitational constant, and angular 
momentum of the black hole, respectively. However, the observations can provide basically 
only two constraints, namely the flux and color temperature (except detailed spectral shapes). 
As a result, the KERRBB solutions are degenerate, so that we need to specify both a and i to 
estimate M and M from the data. Thus, we calculated the best-flt values of M and M, for 
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some representative values of i and a. When we increase a with i fixed, the emergent fiux and 
color temperature both increases, because the last stable orbit gets closer to the black hole. 
Then, the model tries to reduce the color temperature by increasing M, while reducing the fiux 
by lowering M. This a-dependence become more enhanced toward higher inclination, because 
the longitudinal Doppler effects of the disk rotation get stronger. Here, we only considered the 
case of a > (i.e., a prograding accretion disk) and adopted the maximum value of a = 0.998 
(Thorne 1974). 

Acceptable fits have been obtained for a variety of i and a (x^/d.o.f. = 1.07 ~ 1-18); 
the obtained values of M and M are summarised in figure 10, where we also plot the MCD 
model results, in order to illustrate the relativistic effects which are considered in the KERRBB 
model but not in the MCD approximation. In the figure, the Eddington ratio rj is indicated 
with red lines. We regard that the region below the line of r/ = 1 is physically meaningful. 
Thus, the observed X-ray spectrum of Suzaku J1305-4931 can be interpreted as a black hole 
with a ~ 1 and z = 50 -80°, radiating at 77 ^ 1. Especially, a combination of a = 0.998 and i = 60° 
gives the most successful description of the observed spectrum (x^/d.o.f. = 249.3/232), yielding 
M = 3QMq and M = 2.0 x 10^^ g s~^. The fit residuals for these parameters are shown in the 
bottom panel of figure 4. Therefore, we conclude that a nearly extreme Kerr black hole with 
a mass of M = 20 - I3OM0 is one of the most promising candidates for the nature of Suzaku 
J1305-4931. 

Among Galactic/Magellanic black holes, GRS 1915+105 is the only source exhibiting 
both a comparable high Eddington ratio (r/ = 0.2 - 3, Done et al. 2004), and a high spin 
parameter estimated (a~0.7, Middleton et al. 2006). One of the most prominent characteristics 
of GRS 1915+105. especially in the super-Eddington phase (Done et al. 2004), is its rapid 
variability with a timescale from ~ an hour to even down to sub-seconds, presumably due to a 
disk instability (e.g., Greiner et al. 1996; Belloni et al. 1997). Although it is interesting to ask 
whether such a variability is present in Suzaku J1305-4931 as well, we can only say, due to the 
limited statistics (figure 2), that any short-term intensity variation is less than by 10 % and 15 
%, on representative time scale of 1 hour and 10 minutes, respectively. 

We are grateful to all the members of the Suzaku team, for the successful operation 
and calibration. We also thank Prof. R. Narayan for kindly giving us information on the 
KERRBB model. The constructive comments from the anonymous referee greatly helped us 
to improve the quality of the paper. N. I. is supported by the Special Postdoctoral Researchers 
Program in RIKEN. P. G. is a Fellow of the Japan Society for Promotion of Science (JSPS). 
We have made extensive use of the NASA/IPAC Extra galactic Database (NED; the Jet 
Propulsion Laboratory, California Institute of Technology, the National Aeronautics and Space 
Administration) . 
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Table 1. The best-fit spectral parameters of the SRC region for the first observation. 



Parameters Value 
TVh (lO^i cm-2) * 1.6 X 10^ 

hO.26 
-0.06 



kT (keV) 0.34t°-^^ 



Fth (10-13 erg s-i cm'^) t 0.33 

FpL 2.16j;o28 
FpL (10-13 erg s-i cm-2) t 1.6 

xVd.o.f. 43.1/36 

' t^lxed at the Galactic value. 

t Absorption-inclusive 0.5 - 10 keV model flux 

of the thermal component. 

•t Absorption-inclusive 0.5 - 10 kcV model flux 

of the PL component. 
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Table 2. The best-fit parameters for the time-averaged spectrum of Suzaku J1305-4931. 



Parameters 


PL 


MCD 


MCD + PL 


variable-p 


TVh (lO^i cm-2) 


r nr, + 0.43 


r, -. 7+0.27 
^•J^ ' -0.26 


o -14+1.06 


4 19+0.64 

^•^^-0.53 


Tin (keV) 




^•"^-0.05 


1 9r+0.24 

^■^■J-o.ii 


9 01+0.45 

^■>^^-0.23 


Rin (km) * 




77.0t|:? 


lU/ .0_3g 


■5' -"-13.0 


r or p 






1 90+0.51 
"-^ — 1.47 


0.59 ±0.03 


Fx t 


2.4 


2.2 


1.2/1.2 t 


2.2 


xVd.o.f. 


278.7/233 


273.9/233 


235.2/231 


243.5/232 



* innermost disk radius, assuming an inclination ot 60" 

t 0.5 - 10 keV source flux in units of 10~^^ erg cm~^ without correction 
for absorption. 

* flux of PL / MCD components 



Table 3. The spectral variation of Suzaku J1305-4931. 



Parameters 


Period 1 


Period 2 


Period 3 


Period 4 


Exposure (ksec) 


15.8 


13.2 


17.1 


19.9 


Nu (10^1 cm-2) 




2.17 (fix) 




Tin (keV) 


-I 77+0.09 
'-0.08 


^•0^-0.11 


1 RS+OIO 
-L-DO-0.09 


^•"0-0.08 


Rin (km) * 






77.9i^j 


77 7+7.5 
' '-'-6.9 


Fx ^ 


2.9 


1.5 


2.2 


2.0 


Lhoi ^ 


5.9 


3.0 


4.4 


4.1 


xVd.o.L 


135.3/153 


85.7/79 


119.8/132 


143.6/146 



^ Innermost disk radius, assuming an inclination of 60" 
t Absorption-inclusive 0.5 - 10 kcV source flux, in units of 10~^^ erg cm^^ 



■t Bolometric luminosity in units of 10"^^ erg s ^, assuming an inclination of 
60° 
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Fig. 1. XIS FI images of NGC 4945 in 0.5 - 10 keV, obtained in the first (left) and second (right) 
observations. The images are smoothed with a two-dimensional Gaussian of ct = 16" radius. The images 
include background, and are not corrected for vignetting. Regions corresponding to the calibration radio 
active sources at the XIS corners are removed. The overlaid contours show the K-band infrared image 
(Jarrett et al. 2003). The position of the NGC 4945 nucleus is indicated by the filled star, and X-ray sources 
detected with XMM-Newton by black or white crosses. The two circles in the right panel, denoted as SRC 
and BCD, show the integration regions for the source and background signals, respectively, employed to 
study Suzaku J 1305-4931. 



16 



r 


. . . 

Period 1 


1 

T 


Period S ! 


1 

1 Period 3 | 

1 T 


1 M 

1 Period 4 1 : 


r 


1- \ 


1 

+4 


' ' ' ' ! 


1 ' ' ' 1 1 


1 ' ' 1 ' H 


d 

'r 


, , , 


1 
1 


' ! 

+ ,1 


I ... 1 I 


1 + h 

: . . 1 . :n 



50 100 150 

Time [ks] 



Fig. 2. Background-subtracted XIS FI lightcurves of Suzaku J1305-4931, binned into 2 hours (7200 see). 
Tire top and middle panels show the soft (0.5 - 2 keV) and hard (2 - 10 keV) band lightcurves, respectively. 
The hardness of the source, defined as the ratio of the hard band count rate to the soft band one, is shown 
in the bottom panel. The definition of 4 periods, adopted in the study of spectral variations, is indicated 
with dashed lines. 



Source Region (2005/08/22-23) 



MEKAL + PL 




Energy [keV] 



Fig. 3. XIS spectra of the SRC region obtained in the first observation. The FI and BI data are shown 
with the black and red data points, respectively. The best-fit model, consisting of a MEKAL and a PL 
component, are shown with histograms. 



17 



Suzaku J1 305-4931 




12 5 10 

Energy [keV] 



Fig. 4. XIS spectra of Suzaku J1305-4931. The black and red points indicate tlie FI and BI data, respec- 
tively. The best-fit MCD model is shown with the histograms in the first panel, where the contribution 
from the host galaxy is represented by two fixed model components. The second, third, forth and fifth 
panels display the residuals from the PL, MCD, variablc-p disk, and KERRBB (Li ct al. 2005; with i = 60° 
and a = 0.998) fits, respectively. 
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Fig. 5. Spectral energy distribution of Suzaku J1305-4931, in the vF^, form. Only the FI data are shown, 
for clarity. The histograms shows the best-fit MCD model, including the emission from the host galaxy. 
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Fig. 6. The MCD + PL fit to the XIS spectra. Although the same data as in figure 4 are utilized, only 
the FI data are shown for clarity. The blue and green histograms indicate the best-fit MCD and PL 
component, respectively. 
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Fig. 7. Hardness ratio plotted against the total count rate in 0.5 ~ 10 kcV. The blue and green lines 
indicate the prediction for the case of Lboi (i-e., Rm is constant) and Lboi T;^, respectively. The 

source variation is assumed to occur around a reference condition of Ti„ = 1.69 keV and Lboi = 4.4 x 10^^ 
erg s""'^, as determined from the MCD fitting to the time-averaged spectrum shown in figure 4. Each data 
point is an integration for 15 ks. 
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Fig. 8. XIS spectra of Suzaku J1305-4931, obtained in the 4 individual periods defined in figure 2. The 
top-left, top-right, bottom-left, and bottom-right panel displays the spectrum of Period 1, 2, 3 and 4, 
respectively. The histograms in each panel show the best-fit MCD model, including the emission from the 
host galaxy. 
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Fig. 9. The relation between tlie bolometric luminosity iboi and the disk temperature Tin of Suzaku 
J1305-4931 (red) in comparison with those of well-studied ULXs or some luminous cxtragalactic sources 
(blue; Makishima et al. 2000; Kotoku et al. 2000; Mizuno et al. 2001; Miyawaki et al. 2004; Tanaka et al. 
2005; Tsunoda et al. 2006) and galactic black holes (black; Makishima et al. 2000, and reference therein). 
The disk inclination of i = 60° is assumed in the plot. The dotted lines show Lboi expected from the 
standard accretion disk, as iboi oc T-^, assuming Schwarzschild black holes with several representative 
masses. The dashed lines indicate ?y = 1, 0.1 and 0.01. 




Fig. 10. The relation between the mass M and accretion rate M, obtained with the KERRBB model (Li 
et al. 2005) for some representative values of the disk inclination i and spin parameter a. The solid and 
dashed lines indicate constant i and a, respectively. The blue line shows the locus of the MCD model. 
The red lines indicate to r] — 1, 2 and 3. 
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